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t*^ . Dynamical systems theory is especially well-suited for determining the possi- 

ble asymptotic states (at both early and late times) of cosmological models, 
particularly when the governing equations are a finite system of autonomous 
ordinary differential equations. We begin with a brief review of dynamical sys- 
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1 Introduction 

The governing equations of the most commonly studied cosmological models are a 
system of autonomous ordinary differential equations (ODEs). Since our main goal 
is to give a qualitative description of these models, a dynamical systems approach is 
undertaken. Usually, a dimensionless (logarithmic) time variable, r, is introduced so 
that the models are valid for all times (i.e., r assumes all real values). A normalised 
set of variables are then chosen for a number of reasons. First, this normally leads 
to a compact dynamical system. Second, these variables are well-behaved and often 
have a direct physical interpretation. Third, due to a symmetry in the equations, one 
of the equations decouple (in general relativity the expansion is used to normalize the 
variables in ever expanding models whence the Raychaudhuri equation decouples) 
and the resulting simplified reduced system is then studied. The singular points of 
the reduced system then correspond to dynamically evolving self-similar cosmological 
models. More precisely, using the dimensionless time variable and a normalised set 
of variables, the governing ODEs define a flow and the evolution of the cosmological 
models can then be analysed by studying the orbits of this flow in the physical state 
space, which is a subset of Euclidean space. When the state space is compact, each 
orbit will have a non-empty a-limit set and tu-limit set, and hence there will be a 
both a past attractor and a future attractor in the state space. 



1 . 1 Self-similarity 

Self-similar solutions of the Einstein field equations (EFE) play an important role in 
describing the asymptotic properties of more general models. The energy-momentum 
tensor of a perfect fluid given by 

Tab = (/i + P)UaUb + Pgab, (1-1) 

where u"' is the normalized fluid 4-velocity, fi is the density and p is the pressure 
satisfying a linear barotropic equation of state of the form 

p=(7-l)/i, (1.2) 

where 7 is a constant. The existence of a self-similarity of the first kind can be 
invariantly formulated in terms of the existence of a homothetic vector |I|. For a 
general spacetime a proper homothetic vector (HV) is a vector field $, which satisfies 

^^Qfiv = 25'^„, (1.3) 

where g^^^ is the metric and C denotes Lie differentiation along $,. An arbitrary 
constant on the right-hand-side of (1.3) has been rescaled to unity. If this constant is 



^£,guv = 0, then ^ is a Killing vector. A homothetic motion or homothety 
captures the geometric notion of "invariance under scale transformations" . 

Self-similar models are often related to the asymptotic states of more general 
models 0. In particular, self-similar models play an important role in the asymptotic 
properties of spatially homogeneous models, spherically symmetric models, G2 models 
and silent universe models . We will focus on spatially homogeneous models here [^ 
Q We note that the self-similar Bianchi models of relevence below are transitively self- 
similar (in the sense that the orbits of the H4 are the whole spacetime). Self-similar 
spherically symmetric models have been studied by many authors and have been 
recently reviewed in Carr and Coley @]. G2 models, which contain two commuting 
spacelike KV acting orthogonally transitively, have been discussed in Wainwright 
and Ellis (hereafter denoted WE). Other exact homothetic models, including for 
example, plane-symmetric models have been discussed in Kramer et al. |^ and Carr 
and Coley 0. 

Exact solutions 

Let us review some exact self-similar solutions that are of particular importance. 

• Minkowski space (M): 

ds"^ = -dt^ + dx^ + dy"^ + dz^ (1.4) 

i = tdt + x—+y— + z—. (1.5) 

ox ay oz 

In addition to flat Minkowski spacetime, all {k = 0, ±1) Friedmann-Robertson- 
Walker (FRW) models admit a timelike HV in the special case of stiff matter (7 = 2). 
Otherwise, only the k = models admit a HV, and this occurs for all such models in 
which p = (7 — l)/i and hence the scale function has power-law depend on time. 

• fc = 0, FRW (Fo): 

ds"" = -dt"" + t^/^^dx^ + dy"" + dz"") (1.6) 



where 
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p=(7-l)/i; l^=^t h=Qlj- (1-^ 



All transitively self-similar orthogonal spatially homogeneous perfect fluid solu- 
tions (with I < 7 < 2) and spatially homogeneous vacuum solutions are summarized 
in Hsu and Wainwright 0] . In particular, the Kasner vacuum solution is self-similar. 

• Kasner (vacuum) {K): 

ds'^ = -dt^ + t^'PHx^ + t^P^dyS ^ ^2p3^^2 ^^_g) 



with 

Epi = Sp2 = i (1.10) 

f) f) f) 

^ = td, + {l-p,)x— + {l-p2)yg^ + {l-p3)z—.c (1.11) 

1.2 Brief Survey of Techniques in Dynamical Systems 

This section will briefly review some of the results of dynamical systems theory. 
Definition 1 A singular point of a system of autonomous ODEs 

x = f{x) (1.12) 

is a point x G R" such that f{x) =0. 



Definition 2 Let x be a singular point of the DE ( [1.121 ). The point x is called a 



hyperbolic singular point if Re{\i) ^ for all eigenvalues, Aj, of the Jacobian of the 
vector field f{x) evaluated at x. Otherwise the point is called non-hyperbolic. 



Definition 3 Let x(t) = (paif) be a solution of the DE ( |1.12| ) with initial condition 
a;(0) = a. The flow {(7*} is defined in terms of the solution function 0a(^) of the DE 
by 

g^a = (f)a{t). 

Definition 4 The orbit through a, denoted by 7(a) is defined by 

7(a) = {x E R^lx = g^a, for all t E R}. 



Definition 5 Given a DE (|1.12|) in R", a set S* C R" is called an invariant set for 
the DE if for any point a E S the orbit through a lies entirely in S, that is 7(a) C 5*. 

Definition 6 Given a DE ( p..l2| ) in R", with flow {g*}, a subset S C R" is said to 
be a trapping set of the DE if it satisfies: 1. 5* is a closed and bounded set, 2. a E S 
implies that g'^a E S for all t > 0. 

Qualitative analysis of a system begins with the location of singular points. Once 
the singular points of a system of ODEs are obtained, it is of interest to consider the 
dynamics in a local neighbourhood of each of the points. Assuming that the vector 
field f{x) is of class C^ the process of determining the local behaviour is based on 
the linear approximation of the vector field in the local neighbourhood of the singular 
point X. In this neighbourhood 

f{x)^ Df{x){x-x) (1.13) 



where Df{x) is the Jacobian of the vector field at the singular point x. The system 
( |1.13| ) is referred to as the linearization of the DE at the singular point. Each of the 
singular points can then be classified according to the eigenvalues of the Jacobian of 
the linearized vector field at the point. 

The classification then follows from the fact that if the singular point is hyper- 
bolic in nature the fiows of the non- linear system and it's linear approximation are 
topologically equivalent in a neighbourhood of the singular point. This result is given 
in the form of the following theorem: 

Theorem 1: Hartman-Grobman Theorem Consider a DE: x = f{x), where 
the vector field / is of class C^. If a; is a hyperbolic singular point of the DE then 
there exists a neighbourhood of x on which the flow is topologically equivalent to the 
flow of the linearization of the DE at x. 

Given a linear system of ODEs: 

X = Ax, (1.14) 

where A is a matrix with constant coefficients, it is a straightforward matter to 
show that if the eigenvalues of the matrix A are all positive the solutions in the 
neighbourhood oi x = all diverge from that point. This point is then referred to as 
a source. Similarly, if the eigenvalues all have negative real parts all solutions converge 
to the singular point x = 0, and the point is referred to as a sink. Therefore, it follows 
from topological equivalence that if all eigenvalues of the Jacobian of the vector field 
for a non-linear system of ODEs have positive real parts the point is classified as a 
source (and all orbits diverge from the singular point), and if the eigenvalues all have 
negative real parts the point is classified as a sink. 



In most cases the eigenvalues of the linearized system ( |1.13| ) will have eigenvalues 



with both positive, negative and/or zero real parts. In these cases it is important to 
identify which orbits are attracted to the singular point, and which are repelled away 
as the independent variable (usually t) tends to infinity. 

For a linear system of ODEs, ( |1.14| ), the phase space R" is spanned by the eigen- 
vectors of A. These eigenvectors divide the phase space into three distinct subspaces; 
namely: 



^ns) 



The stable subspace E^ = span(si, S2, ...Sr. 
The unstable subspace E^ = span{ui,U2, ■■■Unu) 
The centre subspace E'^ = span(ci,C2, ...Cnc) 

where s, are the eigenvectors who's associated eigenvalues have negative real part, Ui 
those who's eigenvalues have positive real part, and Cj those who's eigenvalues have 
zero eigenvalues. Flows (or orbits) in the stable subspace asymptote in the future to 



the singular point, and those in the unstable subspace asymptote in the past to the 
singular point. 

In the non-linear case, the topological equivalence of flows allows for a similar 
classification of the singular points. The equivalence only applies in directions where 
the eigenvalue has non-zero real parts. In these directions, since the flows are topo- 
logically equivalent, there is a flow tangent to the eigenvectors. The phase space is 
again divided into stable and unstable subspaces (as well as centre subspaces). The 
stable manifold W^ of a singular point is a differential manifold which is tangent to 
the stable subspace of the linearized system [E^). Similarly, the unstable manifold is 
a differential manifold which is tangent to the unstable subspace (-E") at the singular 
point. The centre manifold, W^, is a differential manifold which is tangent to the 
centre subspace E'^. It is important to note, however, that unlike the case of a linear 
system, this centre manifold, W^ will contain all those dynamics not classified by 
linearization (i.e., the non-hyperbolic directions). In particular, this manifold may 
contain regions which are stable, unstable or neutral. The classification of the dynam- 
ics in this manifold can only be determined by utilizing more sophisticated methods, 
such as centre manifold theorems or the theory of normal forms (see |§]). 

Unlike a linear system of ODEs, a non-linear system allows for singular structures 
which are more complicated than that of the singular points, fixed lines or periodic 
orbits. These structures include, though are not limited to, such things as heteroclinic 
and/or homoclinic orbits and non-linear invariant sub-manifolds (for definitions see 
1^). Sets of non-isolated singular points often occur in cosmology and therefore their 
stability will be examined more rigorously. 

Definition 7: A set of non-isolated singular points is said to be normally hy- 
perbolic if the only eigenvalues with zero real parts are those whose corresponding 
eigenvectors are tangent to the set. 

Since by definition any point on a set of non-isolated singular points will have at 
least one eigenvalue which is zero, all points in the set are non-hyperbolic. A set which 
is normally hyperbolic can, however, be completely classified as per it's stability by 
considering the signs of the eigenvalues in the remaining directions (i.e., for a curve, 
in the remaining n — 1 directions) . 

The local dynamics of a singular point may depend on one or more arbitrary 
parameters. When small continuous changes in the parameter result in dramatic 
changes in the dynamics, the singular point is said to undergo a bifurcation. The 
values of the parameter(s) which result in a bifurcation at the singular point can 
often be located by examining the linearized system. Singular point bifurcations will 
only occur if one (or more) of the eigenvalues of the linearized systems are a function 
of the parameter. The bifurcations are located at the parameter values for which the 
real part of an eigenvalue is zero. 



The future and past asymptotic states of a non-linear system may be represented 
by any singular or periodic structure. In the case of a plane system (i.e., in two- 
dimension phase space), the possible asymptotic states can be given explicitly. This 
result is due to the hmited degrees of freedom in the space, and the fact that the flows 
(or orbits) in any dimensional space cannot cross. The result is given in the form of 
the following theorem: 

Theorem 2: Poincare-Bendixon Theorem: Consider the system of ODEs 
X = f{x) on R^, with f ^ C"^, and suppose that there are at most a finite number of 
singular points (i.e., no non-isolated singular points). Then any compact asymptotic 
set is one of the following; 1. a singular point, 2. a periodic orbit, 3. the union of 
singular points and heteroclinic or homoclinic orbits. 

This theorem has a very important consequence in that if the existence of a closed 
(i.e., periodic, heteroclinic or homo-clinic) orbit can be ruled out it follows that all 
asymptotic behaviour is located at a singular point. 

The existence of a closed orbit can be ruled out by many methods, the most 
common is to use a consequence of Green's Theorem, as follows: 

Theorem 3: Dulac's Criterion: If D C i?^ is a simply connected open set 
and V(5/) = gf^(5/i) + ^(5/2) > 0, or (< 0) for all x e D where B is a C^ 
function, then the DE x = f{x) where f E C^ has no periodic (or closed) orbit which 
is contained in D. 

A fundamental criteria of the Poincare-Bendixon theorem is that the phase space is 
two-dimensional. When the phase space is of a higher dimension the requirement that 
orbits cannot cross does not result in such a decisive conclusion. The behaviour in such 
higher-dimensional spaces is known to be highly complicated, with the possibility of 
including such phenomena as recurrence and strange attractors [jlO[ . For that reason. 



the analysis of non-linear systems in spaces of three or more dimensions cannot in 
general progress much further than the local analysis of the singular points (or non- 
isolated singular sets). The one tool which does allow for some progress in the analysis 
of higher dimensional systems is the possible existence of monotonic functions. 

Theorem 4: LaSalle Invariance Principle: Consider a DE x = f{x) on R". 
Let S* be a closed, bounded and positively invariant set of the flow, and let Z he a 
C^ monotonic function. Then for all Xq G S, 

w{xo) c{xE S\Z = 0} 

where w{xo) is the forward asymptotic states for the orbit with initial value xq] i.e., 
a W7-Iimit set [|rT|. 

This principle has been generalized to the following result: 



Theorem 5: Monotonicity Principle (see IT^). Let 0t be a flow on R" with 
5* an invariant set. Let Z : S —>■ H he a C^ function whose range is the interval (a, b), 
where a G RU {— oo}, 6 G RU {c>o} and a < b. If Z is decreasing on orbits in S, then 
for all X eS, 

i-o{x) '^ {s E S \ S\limy^sZ{y) ^ b}, 
<^{x) C {s G S \ S\limy^sZ{y) ^ a}, 

where uj{x) and a{x) are the forward and backward limit set of x, respectively (i.e., 
the w and a limit sets.) 



2 Spatially Homogeneous Perfect Fluid Models 

Many people have studied self-similar spatially homogeneous models, both as exact 
solutions and in the context of qualitative analyses (see WE and Coley |]l3l and ref- 



erences therein). Exact spatially homogeneous solutions were first displayed in early 
papers [|14|; however, it was not until after 1985 that many of them were recognized 
by Wainwright [O and Rosquist and Jantzen |lO, O as being self-similar. The com- 



plete set of self-similar orthogonal spatially homogeneous perfect fluid and vacuum 
solutions were given by Hsu and Wainwright and they have also been reviewed in 
WE. Kantowski-Sachs models were studied by Collins 



Spatially homogeneous models have attracted considerable attention since the gov- 
erning equations reduce to a relatively simple finite-dimensional dynamical system, 
thereby enabling the models to be studied by standard qualitative techniques. Planar 
systems were initially analyzed by Collins ||19|, |2^ and a comprehensive study of gen- 
eral Bianchi models was made by Bogoyavlenski and Novikov [^ and Bogoyavlenski 
p^ and more recently (using automorphism variables and Hamiltonian techniques) 



by Jantzen and Rosquist |T6|, |T^, |2^, |J]. Perhaps the most illuminating approach has 
been that of Wainwright and collaborators [0, ^ |26[ , in which the more physically or 
geometrically natural expansion-normalized (dimensionless) configuration variables 
are used. In this case, the physically admissible states typically lie within a bounded 
region, the dynamical system remains analytic both in the physical region and its 
boundaries, and the asymptotic states typically lie on the boundary represented by 
exact physical solutions rather than having singular behaviour. We note that the 
physically admissible states do not lie in a bounded region for Bianchi models of 
types VIIo, VIII and IX; see WE for details. 

Wainwright utilizes the orthonormal frame method ^ and introduces expansion- 
normalized (commutation function) variables and a new "dimensionless" time vari- 
able to study spatially homogeneous perfect fluid models satisfying p = (7 — l)/i. 



The equations governing the models form an A^-dimensional system of coupled au- 
tonomous ODEs. When the ODEs are written in expansion-normalized variables, 
they admit a symmetry which allows the equation for the time evolution of the ex- 
pansion 6 (the Raychaudhuri equation) to decouple. The reduced A^ — 1-dimensional 
dynamical system is then studied. At all of the singular points of the reduced system, 
9 is proportional to 9"^ and hence all such points correspond to transitively self-similar 
cosmological models . This is why the self-similar models play an important role in 
describing the asymptotic dynamics of the Bianchi models. For orthogonal Bianchi 



models of class A, the resulting reduced state space is five-dimensional ||2^. Orthog- 
onal Bianchi cosmologies of class B were studied by Hewitt and Wainwright |2^ and 
are governed by a five-dimensional system of analytic ODEs with constraints. 



Two perfect-fluid models were studied by Coley and Wainwright |2^. In further 
work, imperfect fluid Bianchi models were studied under the assumption that all 
physical quantities satisfy "dimensionless equations of state" , thereby ensuring that 
the singular points of the resulting reduced dynamical system are represented by 
exact self-similar solutions. Models satisfying the linear Eckart theory of irreversible 
thermodynamics were studied by Burd and Coley [^ and Coley and van den Hoogen 



P^ , those satisfying the truncated causal theory of Israel-Stewart by Coley and van 



den Hoogen [0, and those satisfying the full (i.e., non-truncated) relativistic Israel- 
Stewart theory by Coley et al. [^ . Self-similar solutions also play an important role 
in describing the dynamical behaviour of cosmological models close to the Planck time 
in general relativity with scalar fields [33, ?], in scalar-tensor theories of gravity |P^|, 



and particularly in the low-energy limit in supergravity theories from string theory 
and other higher- dimensional gravity theories. 



2.1 Some Simple Examples 

The expansion scalar 6*, shear scalar a (a^ = \(Tah(^°'^)i and Ricci 3-curvature (orthog- 
onal to u"") ^R, are defined in Ellis [3^. 



Basic equations: 

Raychaudhuri eqn. : 9 = —-9^ — 2cr^ — -{fi + 3p). (2.1) 

o o 

Conservation eqn. : fi = — (/i + p)9. (2.2) 

Generalized Friedmann eqn. : 9'^ = — ^R + 3cr^ + 3 mu. (2.3) 

o 

Ricci identity: for a [35]. (2.4) 

We define the expansion- normalized variables {9 > 0): 

3/i 

il = -rr, density parameter (2.5) 

9"^ 



P = 2v3— , normalized rate of shear (2-6) 

9 

dr 1 

— = -e, [r^-ooast^O+] (2.7) 

A dash (dot) denotes differentiation with respect to T{t). The equation of state is 
given by Eq. (1.2), where 1 < 7 < 2 for normal matter and models with < 7 < | are 
of interest in connection with inflationary models of the universe (see, for example, 
Wald [0). The weak energy condition implies that 



/i>0. 
A. FRW: The metric is given by 



2jr>2 



1 — kr 



+ r'dn 



(2.8) 



where k is the curvature constant. Here a = {9 = 3R/R), and the equations reduce 
to the single ODE: 

n' = {2--f)n{i-n), (2.9) 

where Q > and 1 — ^2 = — | ■^R9~'^, so that f] < 1, = 1, > according to whether 
fc < 0, = 0, > (models open, fiat or closed), respectively. 

• Singular points (7 7^ 2): 

• Q = 1 (-Fq); past attractor (fiat FRW model describes the dynamics 

near the big bang). 

• i7 = (M); future attractor (open models approach Milne form of 

Minlcowski space at later times). 

[For closed models ^ = at the point of maximum expansion, so the models are only 
valid up to this point. The models recollapse and asymptotically approach the fiat 
FRW model Fq to the future] . 

B. Bianchi V: The metric is given by 

ds^ = -dt" + A^{t)dx^ + e^%B\t)dy^ + A\t)B-^{t)dz^), (2.10) 

where 9 = 3^, cr^ = (4 - f), and ^R = -^. The Ricci identity leads to 



& = -9a. (2.11) 

The EFE then reduce to the plane autonomous system: 



a = {2 - 3-f)n{l - Q) + (3^n 
2 



(3' = ^[(37-2)f] + /52-4], (2.12) 
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where 

/3^ + 4fi-4 = --^<0 

(from the generahzed Friedmann eqn.) and 

^ > 0, 

so that phase space is compact. 

• invariant sets: /3 = (FRW, dealt with above). 

P"^ + 4:Q- 4 = (Bianchi I, boundary) 

In the Bianchi I (boundary) case we have a single ODE and the models evolve 
from K to Ft). In addition, the equations above are symmetric about the fi-axis. 

• singular points (fi,/5) [7 7^ 2]: 

• (0, 0) (M) : future attractor 

(1,0) (Fq) : saddle (past attractor for open FRW models, 

intermediate asymptote otherwise). 



(0,2) {K+): 



past attractor (initial Kasner cigar singularity). 



(0,-2) (i^.): 

In the special case 7 = 2, the boundary jS"^ + 4Q — 4 = becomes a line of 
repelling singular points. The corresponding Bianchi I stiff matter models evolve 



from the general self-similar (stiff matter) Bianchi I solution of Jacobs ||38| to the 
Milne (flat) universe [PD[|. 



2.2 Asymptotic States of Bianchi Models 

We now discuss the asymptotic states of Bianchi models, again assuming the linear 
equation of state (1.2). We will summarize the work of Wainwright and Hsu |^ 



and Hewitt and Wainwright [^, who studied the asymptotic states of orthogonal 
spatially homogeneous models in terms of attractors of the associated dynamical 
system for class A and class B models, respectively. Due to the existence of monotone 
functions, it is known that there are no periodic or recurrent orbits in class A models. 
Although "typical" results can be proved in a number of Bianchi type B cases, these 
are not "generic" due to the lack of knowledge of appropriate monotone functions. 
In particular, there are no sources or sinks in the Bianchi invariant sets B^ (VIII) or 
5±(IX). 

The key results are as follows: 
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• A large class of orthogonal spatially homogeneous models (including all class 
B models) are asymptotically self-similar at the initial singularity and are ap- 
proximated by exact perfect fluid or vacuum self-similar power-law models. 
Examples include self-similar Kasner vacuum models or self-similar locally ro- 
tationally symmetric (class III) Bianchi type II perfect fluid models |19|, 0, p9| . 

However, this behaviour is not generic; general orthogonal models of Bianchi types 
IX and VIII have an oscillatory behaviour with chaotic-like characteristics, with the 
matter density becoming dynamically negligible as one follows the evolution into the 
past towards the initial singularity. Ma and Wainwright |^ show that the orbits of 
the associated cosmological dynamical system are negatively asymptotic to a lower 
two-dimensional attractor. This is the union of three ellipsoids in R^ consisting of 
the Kasner ring joined by Taub separatrices; the orbits spend most of the time near 
the Kasner vacuum singular points. Clearly the self-similar Kasner models play a 
primary role in the asymptotic behaviour of these models. 

• Exact self-similar power-law models can also approximate general Bianchi mod- 
els at intermediate stages of their evolution (e.g., radiation Bianchi Vllh models 
^). Of special interest are those models which can be approximated by an 
isotropic solution at an intermediate stage of their evolution (e.g., those models 
whose orbits spend a period of time near to a fiat Friedmann singular point). 

This last point is of particular importance in relating Bianchi models to the real 
Universe, and is discussed further in general terms in WE (see, especially. Chapter 
15) and specifically for Bianchi Vllh models in Wainwright et al. |^. In particular, 
the fiat Friedmann singular point is universal in that it is contained in the state 
space of each Bianchi type. Isotropic intermediate behaviour has also been found in 
tilted Bianchi V models |42|, and it appears that many tilted models have isotropic 



intermediate behaviour (see WE). 

• Self-similar solutions can describe the behaviour of Bianchi models at late times 
(i.e., as t ^ oo). Examples include self-similar fiat space and self-similar homo- 



geneous vacuum plane waves [19, 15 . 



All models expand indefinitely except for the Bianchi type IX models. The question 
of which Bianchi models can isotropize was addressed in the famous paper by Collins 
and Hawking [^, in which it was shown that, for physically reasonable matter, the 



set of homogeneous initial data that give rise to models that isotropize asymptotically 
to the future is of zero measure in the space of all homogeneous initial data (see also 
WE). 
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All vacuum models of Bianclii (B) types IV, V, Vlh and Vllh are asymptotic to plane 

wave states to the future. Type V models tend to the Milne form of fiat spacetime 

26| . Typically, and perhaps generically |26], non- vacuum models are asymptotic in 



the future to either plane-wave vacuum solutions [^ or non-vacuum Collins type Vlh 
solutions |T9| . 

Bianchi (A) models of types VIIo (non-vacuum) and VIII expand indefinitely but 
are found to have oscillatory (though non-chaotic) behaviour in the Weyl curvature 
(see, for example, [0). Bianchi type IX models obey the "closed universe recoUapse" 
conjecture p5|. All orbits in the Bianchi invariant sets B^{VIIo) (fi > 0), B^{VIII) 



and B^{IX) are positively departing; in order to analyse the future asymptotic states 
of such models it is necessary to compactify phase-space. The description of these 
models in terms of conventional expansion-normalized variables is only valid up to 
the point of maximum expansion (where ^ = 0). 



3 Scalar Field Models 



3.1 Background 



A variety of theories of fundamental physics predict the existence of scalar fields | J6, 
47| , ^ , motivating the study of the dynamical properties of scalar fields in cosmology. 
Indeed, scalar field cosmological models are of great importance in the study of the 



early universe, particularly in the investigation of inflation ^\ ^, ^. Recently 
there has also been great interest in the late-time evolution of scalar field models. 
'Quintessential' scalar field models (or slowly decaying cosmological constant models) 
pT| , |52[ give rise to a residual scalar field which contributes to the present energy- 
density of the universe that may alleviate the dark matter problem and can predict 
an effective cosmological constant which is consistent with observations of the present 
accelerated cosmic expansion (ESl Q]. 



Models with a self-interaction potential with an exponential dependence on the 
scalar field, 0, of the form 

V = Ae^*, (3.1) 

where A and k are positive constants, have been the subject of much interest and 

arise naturally from theories of gravity such as scalar-tensor theories or string theories 

W^ . Recently, it has been argued that a scalar field with an exponential potential 



is a strong candidate for dark matter in spiral galaxies |^ and is consistent with 



observations of current accelerated expansion of the universe [^ 



A number of authors have studied scalar field cosmological models with an ex- 
ponential potential within general relativity. Homogeneous and isotropic FRW mod- 
els were studied by Halliwell [^ using phase-plane methods. Homogeneous but 
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anisotropic models of Bianchi types I and III (and Kantowski-Sachs models) were 
studied by Burd and Barrow |^, Bianchi type I models were studied by Lidsey |^ 



and Aguirregabiria et al. [Q , and Bianchi models of types III and VI were studied 
by Feinstein and Ibanez [^. A qualitative analysis of Bianchi models with k"^ < 2 
(including standard matter satisfying standard energy conditions) was completed by 
Kitada and Maeda |6^ . The governing differential equations in spatially homogeneous 



Bianchi cosmologies containing a scalar field with an exponential potential reduce to 
a dynamical system when appropriate expansion- normalized variables are defined. 
This dynamical system was studied in detail in |B^ (where matter terms were not 
considered) . 

One particular solution that is of great interest is the flat, isotropic power-law 
inflationary solution which occurs for k"^ < 2. This power-law inflationary solution 
is known to be an attractor for all initially expanding Bianchi models (except a 



subclass of the Bianchi type IX models which will recollapse) |^2|, ^. Therefore, all 
of these models inflate forever; there is no exit from inflation and no conditions for 
conventional reheating. 

Recently cosmological models which contain both a scalar field with an exponential 
potential and a barotropic perfect fluid with a linear equation of state given by (1.2), 
where 7 is in the physically relevant range 2/3 < 7 < 2, have come under heavy 
analysis. One class of exact solutions found for these models has the property that 
the energy density due to the scalar field is proportional to the energy density of the 
perfect fluid, and hence these models have been labelled matter scaling cosmologies 
0, ^, ^. These matter scaling solutions are spatially flat isotropic models and are 
known to be late-time attractors (i.e., stable) in the subclass of flat isotropic models 
|65| and are clearly of physical interest. In addition to the matter scaling solutions. 



curvature scaling solutions [^] and anisotropic scahng solutions [^ are also possible. 
A comprehensive analysis of spatially homogeneous models with a perfect fluid and 
a scalar field with an exponential potential has recently been undertaken [^ . 

Although the exponential potential models are interesting models for a variety of 
reasons, they have some shortcomings as inflationary models |^, ^ . While Bianchi 
models generically asymptote towards the power-law inflationary model in which 
the matter terms are driven to zero for /c^ < 2, there is no graceful exit from this 
inflationary phase. Furthermore, the scalar field cannot oscillate and so reheating 



cannot occur by the conventional scenario. In recent work |]69| interaction terms were 
included, through which the energy of the scalar field is transferred to the matter 
fields. These terms were found to affect the qualitative behaviour of these models 
and, in particular, lead to interesting inflationary behaviour. 
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3.2 Isotropisation 

In the famous paper by Collins and Hawking |Q it was proven that within the 
set of spatially homogeneous cosmological models (which satisfy reasonable energy 
conditions) those which approach isotropy at infinite times is of measure zero; that is, 
in general anisotropic models do not isotropize as they evolve to the future. Since we 
presently observe the universe to be highly isotropic, we therefore need an explanation 
of why our universe has evolved the way it has. This problem, known as the isotropy 
problem, can be easily solved with an idea popularized by Guth ||5^. If the early 



universe experiences a period of inflation, then all anisotropics are essentially pushed 
out of our present observable light-cone and are therefore not presently observed. The 
Cosmic No-Hair Conjecture asserts that under appropriate conditions, any universe 
model will undergo a period of inflation and will consequently isotropize. 

A signifigant amount of work on the Cosmic No-hair Conjecture has already been 
done for spatially homogeneous (Bianchi) cosmologies |3^, [7^, |7T|, |7^, |7^. For in- 
stance, Wald [^ has proven a version of the Cosmic No-Hair Conjecture for spatially 



homogeneous spacetimes with a positive cosmological constant; namely, he has shown 
that all initially expanding Bianchi models asymptotically approach a spatially ho- 
mogenous and isotropic model, except the subclass of Bianchi type IX models which 
recoUapse. 

Anisotropic models with scalar fields and with particular forms for the scalar field 
potential have also been investigated. Heusler [^ has analyzed the case in which 
the potential function passes through the origin and is concave up and, like Collins 
and Hawking |53] has found that the only models that can possibly isotropize to the 



future are those of Bianchi types I, V, VII and IX. 

As noted above, Kitada and Maeda [^, Q have proven that li k < a/2, then all 
initially expanding Bianchi models except possibly those of type IX must isotropize. 
Let us consider what happens in the case k > \^. In Ibaiiez et al. |^ it was 



proven, using results from Heusler's paper |7ll], that the only models that can possibly 



isotropize when k > ^/2 are those of Bianchi types I, V, VII, or IX. Since the Bianchi 
I, V and VHq models are restricted classes of models, the only general spatially 
homogeneous models that can possibly isotropize are consequently of types Vllh or 
IX. Here we shall study the possible isotropization of the Bianchi type Vllh models 
when k > \f2. 



3.2.1 The Bianchi VII,, Equations 

The Bianchi type VII/j models belong to the Bianchi type B models as classified by 



Ellis and MacCallum [ffl. Hewitt and Wainwright \lt\ have derived the equations 



describing the evolution of the general Bianchi type B models. We shall utilize these 
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equations, adjusted so that they describe a model with a minimally coupled scalar 
field (f) with an exponential potential 1^(0) = Ke'^'^. The energy-momentum tensor 
describing a minimally coupled scalar field is given by 

Tab = 4>;a(p;b " 9ab (^<P;c4>''' + V{(f))j , 

where, for a homogeneous scalar field, = (f){t). In this case we can formally treat the 



energy-momentum tensor as a perfect fluid with velocity vector u"" = (p'""/ J— 
where the energy density, /x^, and the pressure, p^, are given by 

1 • 



1 



/i</, = -r + Vi(l)), (3.2) 



202-\/(0). (3.3) 



Our variables are the same as those used by Hewitt and Wainwright ||2^, with the 
addition of 

to describe the scalar field. We note that Q^ = Sfi^/O"^ = x^ + y"^. 
The dimensionless evolution equations are then p6| : 

SV = (g-2)S+-2iV, (3.5) 

(3.6) 
(3.7) 
(3.8) 
(3.9) 

(3.10) 

where the prime denotes differentiation with respect to the new dimensionless time 



(3.11) 
y\ (3.12) 

(3.13) 

iV > 0. (3.14) 



S' 


= 2(g - 2)S - 4AA^+ - 4S+A, 


A' 


= 2(g + S+ - 1)A + 2(S - iV)Ar+, 


A 


= 2(g + 2S+)i, 


% 


= (g + 2S+)Ar+ + 6A, 

1 — 


x' 


= {q-2)x-Jhy\ 



r, where 4^ = §, and 








N = 


liNl-iA), 




Q = 


2S^ + 2S + 2x^ 


There also exists the constraint. 






EiV 


'-a^-aj:1 = 


and the equations are 


subject to the conditions 




i>0, 


S>0, 
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The generalized Friedinann equation, written in diniensionless variables, becomes 

y'^ = l-T.l-t-A-N-x^, (3.15) 

which serves to define j/, and the evolution of y is governed by 

y' = {q + l + \hx)y. (3.16) 



Eqs. ( |3.1CI| ) and ( |3.16| ) are equivalent to the Klein-Gordon equation (f) + 9(f) + V\(j)) = 



written in dimensionless variables. The parameter i = ji defines the group parameter 
h in the Bianchi Vllh models. 

The variables S_|_ and E describe the shear anisotropy. The variables A, N^ and 
N describe the spatial curvature of the models. The variable A describes the relative 
orientation of the shear and spatial curvature eigenframes. 

We are not interested here in the complete qualitative behaviour of the cosmolog- 
ical models [^ but simply whether the Bianchi Vllh models isotropize to the future 
when k > \/2. This question can be easily answered by examining the stability of the 
isotropic singular points of the above six-dimensional dynamical system ( p.5| - p.lO| ). 



3.2.2 Stability Analysis 

All of the isotropic singular points lie in the invariant set FRW defined by {S+ = 
0, S = 0, A = 0, A = 0}. Therefore, we shall find all of the isotropic singular points 



and determine whether any are stable attractors or sinks |75 
The singular point 

S+ = 0,I] = 0,A = 0,v4 = l,A+ = v^,s = (3.17) 

implies y = and represents the negatively curved Milne vacuum model. The lin- 
earization of the dynamical system in the neighborhood of this singular point has 
eigenvalues 

0, 2, -2, -4, -2 + 4v^, -2 + 4v^. (3.18) 

Therefore, this singular point is a saddle. 
The singular point (s) 

S+ = 0, S = 0, A = 0, A = 0, A+ = 0, X = 1 or - 1 (3.19) 

imply that y = and represent fiat non-infiationary FRW model(s). The eigenvalues 
in both cases are 

0,0,2,2,4,6 + V6A;. (3.20) 
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These singular points are unstable with an unstable manifold of at least dimension 4. 
The singular point 



0, S = 0, A = 0, A = 0, A^i 



0,a; 



6 ^ 



implies that y 



1 



fc2 



. The eigenvalues are 



i_^, 1_^, A;2 - 6, A;^ - 4, A;2 - 2, ^-^ 



(3.21) 



(3.22) 



For k < v2, this singular point represents the usual power-law inflationary attractor. 
If \/2 < k < 2, then the singular point has an unstable manifold of dimension 4. If 
2 < k < v^, then the singular point has an unstable manifold of dimension 3. This 
singular point does not exist if A; > y^- 



The singular point 



0,S = 0,A = 0,A 



F' 



N^ 



2 



1 - — l,a; = - 



3k 



(3.23) 



denoted F, implies that y = ^^ and represents a non- inflationary negatively curved 
FRW model. The eigenvalues are 



-2 + ^^{k^-A£k^ + 8£)+E, 
-2 - ^^{k^-A£k^ + M)+E, 

rv 



-2 + ^^{k^ - 4£P + M) - E,(3.24) 
V2„„ .„,.„„ ^_ 



-2-^J{k^-A£k^ 
k ^ 



where 



E = J{k^ + Aik^ - Uf + 32£(2 - A;2 



After some algebra it can be shown that iik > a/2 (note that £ > in the Bianchi VII/i 
models) then all of the eigenvalues have negative real parts. Therefore, if /c > y^, 
then this singular point is a stable attractor. (Note that this singular point does not 
exist if A; < \/2.) In other words, there exists an open set of initial conditions in the 
set of anisotropic Bianchi VII/i (with a scalar fleld and exponential potential) initial 
data for which the corresponding cosmological models asymptotically approach an 
isotropic and negatively curved FRW model. 
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3.2.3 Discussion 

We have shown that within the set of all spatially homogeneous initial data, there 
exists an open set of initial data describing the Bianchi type VII^ models (having a 
scalar field with an exponential potential and k > a/2) such that the models approach 
isotropy at infinite times. This compliments the results of Kitada and Maeda [^, |62| , 
who showed that all ever-expanding spatially homogeneous models (including the 
Bianchi Vllh models) with k < \pl approach isotropy to the future. In other words, 
there exists a set of spatially homogeneous initial data of non-zero measure for which 
models will isotropize to the future for all positive values of k. Of course, there 
also exists a set of spatially homogeneous initial data of non-zero measure for which 
models will not isotropize to the future when k > \pl (e.g., the Bianchi VIII models). 

\ik < -\/2, then all models will inflate as they approach the power-law inflationary 
attractor represented by Eq. (|3.21| ). For k > -\/2, the stable singular point F, 



given by Eq. ( |3.23| ), which does not exist for k < -\/2, is isotropic and resides on 
the surface q = 0. This means that the corresponding exact solution is marginally 
non- inflationary. However, this does not mean that the corresponding cosmological 
models are not inflating as they asymptotically approach this singular state. As 
orbits approach F they may have g < or g > (or even q = 0) and consequently 
the models may or may not be inflating. If they are inflating, then the rate of inflation 



is decreasing as F is approached (i.e., g — > 0). When \/2 < k < J 8/ 3, we find that 
F is node-like, hence there is an open set of models that inflate as they approach F 
and an open set which do not. When k > y 8/3, F is found to be spiral-like, and so 
it is expected that orbits experience regions of both g < and g > as they wind 
their way towards F. As in Kitada and Maeda |7^ , |62| , the inclusion of matter in the 
form of a perfect fluid is not expected to change the results of the analysis provided 
the matter satisfies appropriate energy conditions. 



3.3 Stability of Matter Scaling Solutions 

Spatially homogeneous scalar field cosmological models with an exponential potential 
and with barotropic matter may also be important even if the exponential potential 
is too steep to drive inflation. For example, there exist 'scaling solutions' in which the 
scalar field energy density tracks that of the perfect fluid (so that at late times neither 
field is negligible) |B^. In particular, in |^ a phase-plane analysis of the spatially flat 



FRW models showed that these scaling solutions are the unique late-time attractors 
whenever they exist. The cosmological consequences of these scaling models have 
been further studied in [Q. For example, in such models a significant fraction of the 
current energy density of the Universe may be contained in the homogeneous scalar 
field whose dynamical effects mimic cold dark matter; the tightest constraint on these 



cosinological models comes from primordial nucleosynthesis bounds on any such relic 



density 



Clearly these matter scaling models are of potential cosmological significance. It 
is consequently of prime importance to determine the genericity of such models by 
studying their stability in the context of more general spatially homogeneous models. 



3.3.1 The Matter Scaling Solution 



The governing equations for a scalar field with an exponential potential V = V^e^'^ 
evolving in a flat FRW model containing a separately conserved perfect fluid which 
satisfies the barotropic equation of state 



P7 = (7 - l)/^7' 



where < 7 < 2 here, are given by 

e = 

= 

subject to the Friedmann constraint 



"2(7/^7 + 
-66 - kV, 



^2 = 3(/x, + ^02^\/), 



(3.25) 

(3.26) 
(3.27) 



(3.28) 



where an overdot denotes ordinary differentiation with respect to time t, and units 
have been chosen so that SvrG = 1. We note that the total energy density of the 
scalar field is given by Eq. 



Defining x and y by Eq. ( p.4|) and again using the logarithmic time variable, r, 
Eqs. ( p. 251 ) - ( |3.27| ) can be written as the plane-autonomous system [p5[] : 



X 



—3a; 



-y 

T 



'^V + ^x[2a;2 + 7(l-x2-y')], 



kx^2x^ + 7(1 



X 



(3.29) 
(3.30) 



where 



fi 



3/i^ 



VIa. = 



3/i0 



x^ + 2/^; fi + l^ri 



^2 ' '^- Q2 
which implies that < a;^ + y^ < 1 for fi > 0, so that the phase-space is bounded 



(3.31) 



A qualitative analysis of this plane-autonomous system was given in |B^ . The well- 
known power-law infiationary solution for k'^ < 2 M7\, p2| corresponds to the singular 
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point X = -k/VQ, y = (l- A;V6)i/2 (^^ = 1, fi = 0) of the system (^1291) /(^30|), 
which is shown to be stable (i.e., attracting) for k'^ < 87 in the presence of a barotropic 
fluid. Previous analysis had shown that when k^ < 2 this power-law inflationary 
solution is a global attractor in spatially homogeneous models in the absence of a 
perfect fluid (except for a subclass of Bianchi type IX models which recollapse). 

In addition, for 7 > there exists a scaling solution corresponding to the singular 
point 

x = xo = -Jl^, y = yo= [3(2 -7)7/2^2]^, (3.32) 



whenever k"^ > 87. The linearization of system ( |8.29| )/ (|3.8CI| ) about the singular point 



( p.32| ) yields the two eigenvalues with negative real parts 



- ^ (2 - 7) ± ^V(2 - 7)[24y - P(97 - 2)] (3.33) 

when 7 < 2. The singular point is consequently stable (a spiral for k'^ > 247^/(97 — 2), 
else a node) so that the corresponding cosmological solution is a late-time attractor 
in the class of flat FRW models in which neither the scalar field nor the perfect fluid 
dominates the evolution. The effective equation of state for the scalar field is given 
by 

'0 — 2 I 2 7) 

/^0 xo + yo 

which is the same as the equation of state parameter for the perfect fluid. The 
solution is referred to as a matter scaling solution since the energy density of the 
scalar field remains proportional to that of the barotropic perfect fluid according to 
fi/fi(f) = fc^/37 — 1 [^. Since the scaling solution corresponds to a singular point of 



the system ( p.29| )/( 3.3C ) we note that it is a self-similar cosmological model 



3.3.2 Stability of the Matter Scaling Solution 

Let us study the stability of the matter scahng solution with respect to anisotropic 



and curvature perturbations within the class of spatially homogeneous models ||76|, p2| . 

Bianchi I models 

In order to study the stability of the scaling solution with respect to shear pertur- 
bations we shall first investigate the class of anisotropic Bianchi I models, which are 
the simplest spatially homogeneous generalizations of the flat FRW models and have 
non-zero shear but zero three-curvature. The governing equations in the Bianchi I 
models are Eqs. ( p.26| ) and ( |3.27|) , and Eq. ( p.28|) becomes 



9'' = 3[fi^ + l-<P^ + V)+J:\ (3.34) 
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where T? = SSg-R ^ is the contribution due to the shear, where Sq is a constant and 
R is the scale factor. Eq. (|3.25|) is replaced by the time derivative of Eq. (|3.34| ). 

Using the definitions ( p.4|) and ( p.31j ) we can deduce the governing ODEs. Due to 
the E^ term in ( |3.34|) we can no longer use this equation to substitute for /i^ in the 
remaining equations, and we consequently obtain the three-dimensional autonomous 
system: 



x' -- 


= -3x - ^hy' + ^x[2 + (7 - 2)Q 


y' - 


= lyUlkx + 2 + i^-2)Q-2y' 


n' -- 


= 3n{{^-2){n-l)-2y'}, 



2y' 



where Eq. ( |3.34|) yields 



1-n-x^ -y' = E'e- 

so that we again have a bounded phase-space. 



>0, 



(3.35) 

(3.36) 
(3.37) 

(3.38) 



The matter scaling solution, corresponding to the flat FRW solution, is now rep- 
resented by the singular point 



X 



xo, y = yo, 



Q 



A;2 



(3.39) 



The linearization of system ( |3.35|) - ( p.37| ) about the singular point (|3.39| ) yields three 
eigenvalues, two of which are given by (|3.33|) and the third has the value —3(2 — 7), 
all with negative real parts when 7 < 2. Consequently the scaling solution is stable 
to Bianchi type I shear perturbations. 

Curved FRW models 

In order to study the stability of the scahng solution with respect to curvature 
perturbations we shall first consider the class of FRW models, which have curvature 
but no shear. Again Eqs. ( p.26|) and (|3.27|) are valid, but in this case Eq. ( p.28|) 
becomes 

9^ = 3(/i^ + i^' + V) + K, (3.40) 

where K = —9kR~^ and A; is a constant that can be scaled to 0, ±1. Eq. ( |3.25| ) is 
again replaced by the time derivative of Eq. ( p.40| ). 



As in the previous case we cannot use Eq. ( p.40| ) to replace /i^, and using the 
definitions (|3.4|) and ( p.31|) we obtain the three-dimensional autonomous system: 



X 



-3x 



3, 2 3 
2% +2^ 



7-^)^ + ^(1 + 2.^ 



(3.41) 
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V' = ^y{\l^kx+(j~^^Q+^{l + 2x'-y')}, (3.42) 




n' = 3n\{^-^]in-i) + '^i2x'-y')\, (3.43) 



where 



1-Q-x'^ -y^ = KH-^. (3.44) 

The phase-space is bounded for i^ = or i^ < 0, but not for i^ > 0. 



The matter scahng solution again corresponds to the singular point ( p.39|) . The 



linearization of system ( |3.41| ) - ( p.43| ) about this singular point yields the two eigen- 



values with negative real parts given by ( p.33| ) and the eigenvalue (37 — 2). Hence 
the scaling solution is only stable for 7 < |. For 7 > | the singular point ( ^.39| ) 



is a saddle with a two-dimensional stable manifold and a one-dimensional unstable 
manifold. 

Consequently the scaling solution is unstable to curvature perturbations in the 
case of realistic matter (7 > 1); i.e., the scaling solution is no longer a late-time 
attractor in this case. However, the scaling solution does correspond to a singular 
point of the governing autonomous system of ODEs and hence there are cosmological 
models that can spend an arbitrarily long time 'close' to this solution. Moreover, 
since the curvature of the Universe is presently constrained to be small by cosmolog- 
ical observations, it is possible that the scaling solution could be important in the 
description of our actual Universe. That is, not enough time has yet elapsed for the 
curvature instability to have effected an appreciable deviation from the flat FRW 
model (as in the case of the standard perfect fluid FRW model). Hence the scaling 
solution may still be of physical interest. 

Bianchi Vllh models 

To further study the significance of the scaling solution it is important to deter- 
mine its stability within a general class of spatially homogeneous models such as the 
(general) class of Bianchi type Vllh models (which are perhaps the most physically rel- 
evant models since they can be regarded as generalizations of the negative-curvature 
FRW models). The Bianchi Vllh models are sufficiently complicated that a simple 
coordinate approach (similar to that given above) is not desirable. In subsection 3.2.1 
the Bianchi Vllh spatially homogeneous models with a minimally coupled scalar field 
with an exponential potential (but without a barotropic perfect fluid) were studied by 
employing a group-invariant orthonormal frame approach with expansion-normalized 
state variables governed by a set of dimensionless evolution equations (constituting 
a 'reduced' dynamical system) with respect to a dimensionless time subject to a 
non-linear constraint [Q. A barotropic perfect fluid can easily be included [ p^ . 

The reduced dynamical system is seven- dimensional (subject to a constraint). 
The scaling solution is again a singular point of this seven-dimensional system. This 
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singular point, which only exists for k'^ > 87, has two eigenvalues given by ( p.33|) 



which have negative real parts for 7 < 2, two eigenvalues (corresponding to the 
shear modes) proportional to (7 — 2) which are also negative for 7 < 2, and two 
eigenvalues (essentially corresponding to curvature modes) proportional to (87 — 2) 
which are negative for 7 < § and positive for 7 > | |]T6[. The remaining eigenvalue 
(which also corresponds to a curvature mode) is equal to 87 — 4. Hence for 7 < | 
{k"^ > 87) the scaling solution is again stable. However, for realistic matter (7 > 1) 
the corresponding singular point is a saddle with a four- or five-dimensional stable 
manifold (depending upon whether 7 > 4/8 or 7 < 4/8, respectively). 



4 String Models 



There has been considerable interest recently in the cosmological implications of string 
theory. String theory introduces significant modifications to the standard, hot big 
bang model based on conventional Einstein gravity. Early-universe cosmology pro- 
vides one of the few environments where the predictions of the theory can be quanti- 
tatively investigated. The evolution of the very early universe below the string scale is 
determined by ten-dimensional supergravity theories ||46| , |77| . All theories of this type 
contain a dilaton, a graviton and a two-form potential in the Neveu-Schwarz/Neveu- 
Schwarz (NS-NS) bosonic sector. If one considers a Kaluza-Klein compactification 
from ten dimensions onto an isotropic six-torus of radius e^, the effective action is 
given by 



S 



d x\ 



-ge 



R + (V*)' - 6 (V/3)' - ^H,„.,H^''^ 



Ifii/X-^ 



(4.1) 



where the moduli fields arising from the compactification of the form-fields on the 
internal dimensions and the graviphotons originating from the compactification of 
the metric have been neglected [Q. In Eq. (|4.1| ), R is the Ricci curvature of the 
spacetime with metric g^^ and g = detg^^, the dilaton field, $, parametrizes the string 
coupling, (yf^ = e*, and H^^x = d[^iBy\] is the field strength of the two-form potential 
Bav The volume of the internal dimensions is parametrized by the modulus field, j3. 



In four dimensions, the three-form field strength is dual to a one-form: 



(4.2) 



where e^'^^'^ is the covariantly constant four-form. In this dual formulation, the field 
equations (FE) can be derived from the action 



S = d*Xy/—ge * 



R + {V<l>f - 6 (V/3)' - ^e^* {Vaf 



(4.8) 



where a is interpreted as a pseudo-scalar 'axion' field [179 1. 
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It can be shown that the action ( [4.3|) is invariant under a global SL(2, R) trans- 
formation on the dilaton and axion fields [^. The general FRW cosmologies derived 
from Eq. ( |4.3|) have been found by employing this symmetry [^]. However, the sym- 
metry is broken when a cosmological constant is present ^l|] and the general solution 
is not known in this case. The purpose here is to determine the general structure 
of the phase space of solutions for the wide class of string cosmologies that contain 
a cosmological constant in the effective action. This is particularly relevant in light 
of recent high redshift observations that indicate a vacuum energy density may be 
dominating the large-scale dynamics of the universe at the present epoch |]52[ . 



A cosmological constant may arise in a number of different contexts and we con- 
sider a general action of the form 



S 



d^Xy 



-9 



,-* 



R + iv^f -eivpy 



1 

-e 
2 



2$ 



(Vaf - 2A 



A 



R 



(4.4) 



The constant, A, is determined by the central charge deficit of the string theory and 
may be viewed as a cosmological constant in the gravitational sector of the theory. 
In principle, it may take arbitrary values if the string is coupled to an appropriate 
conformal field theory. Such a term may also have an origin in terms of the reduction 
of higher degree form-fields |Q. The constant, Ar, represents a phenomenological 
cosmological constant in the matter sector. Since it does not couple directly to the 
dilaton field, it may be viewed in a stringy context as a Ramond-Ramond (RR) degree 
of freedom (a 0-form) [|^. Such a cosmological constant may also be interpreted as 
the potential energy of a scalar field that is held in a false vacuum state. 

We shall include the combined effects of the axion, modulus and dilaton fields, 
thereby extending previous qualitative analyses where one or more of these terms was 
neglected |]85| , ^ |87| , |88| . A full stability analysis can be performed for all models by 



rewriting the FE in terms of a set of compactified variables. As usual, units in which 
c = SttG = 1 will be utihzed throughout. 



4.1 Cosmological Field Equations 



When Ar = 0, the spatially fiat FRW cosmological FE derived from action (Hl:.4|) are 
given by 

2a - 2d0 - (T'e2^+^" : 
2y3 - 02 - 3a^ - 6$' + -a^e^^^^^ + 2A : 

a + a {if + 6d) -- 

where yj = $ — 3a defines the 'shifted' dilaton field, R{t) = e" is the scale factor of 
the universe and a dot denotes differentiation with respect to cosmic time, t. The 






(4.5) 





(4.6) 





(4.7) 


0, 


(4.8) 
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generalized Friedinann constraint equation is 



Sd" - 0' + 6/3" + ^d2e2'^+^° + 2A = 0. (4.9) 



A number of exact solutions to Eqs. (|4.5| )- (|4.9| ) are known when one or more of the 
degrees of freedom are trivial; these solutions lie in the invariant sets of the full phase 
space. The 'dilaton-vacuum' solutions, where only the dilaton field is dynamically 
important, are given in [^; there is a curvature singularity in these solutions at 
t = 0. In the pre-big bang inflationary scenario, the pre-big bang phase corresponds 
to the range t < and the post-big bang phase to the solution for t > 0. The 
'dilaton-moduli- vacuum' solutions have d = A = 0. The general solution with A = 



is the 'dilaton-moduli-axion' solution [80|; this cosmology asymptotically approaches 



a dilaton-moduli-vacuum solution in the limits of high and low spacetime curvature. 
The axion field induces a smooth transition between these two power-law solutions. 
The solutions where only the axion field is trivial and A > are specific cases of the 
'rolling radii' solutions found by Mueller |9^. The corresponding solutions for A < 



are related by a redefinition. Finally, there exists the 'linear dilaton-vacuum' solution 
where A > [P^. This solution is static and the dilaton evolves linearly with time. 



4.2 Qualitative Analysis of the NS-NS Fields 

For an arbitrary central charge deficit, the FE (|4.5| )- (|4.9| ) may be written as an 
autonomous system of ODEs: 

(4.10) 
(4.11) 
(4.12) 
(4.13) 

(4.14) 



A^ = 6/5', p = a'e"^+''", ^ = 0, h = a, (4.15) 

and h = 6/3, where 6 is the expansion scalar defined earlier. The variable p may 
be interpreted as the effective energy density of the pseudo-scalar axion field [p8| . 
It follows from Eq. ( [4.11|) that ip is a monotonically increasing function of time and 



h 


= ip'^ + hi!- 


-3h'^-N 


-2A 


^ 


= 3h^ + N 






N 


= 2Ni) 






P 


= —Qhp 






3h' 


- %l)'^ + N + 


ip + 2A: 


= 0, 


e new variables 






5 


p = ^26^^+^", 


^ = 


0, 



this implies that the singular points of the system of ODEs must be located either at 
zero or infinite values oi ip- In addition, due to the existence of a monotone function, 
it follows that there are no periodic or recurrent orbits in the corresponding phase 
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space [0, |9^. The sets A = and p = are invariant sets. In particular, the exact 
solution for A = divides the phase space and the orbits do not cross from positive 
to negative A. 

We must consider the cases where A < and A > separately. In the case 
where the central charge deficit is negative, A < 0, it proves convenient to employ 
the generalized Friedmann constraint equation ([4.14|) to eliminate the modulus field. 



We may compactify the phase space by normalizing with \/i\}^ — 2A and we define a 

new time variable by 

d _ \ d 

dT ^ y^i;-' -2A Jt' ^^'^^^ 

The governing equations reduce to a three-dimensional system of autonomous ODEs. 
The singular points all lie on one of the two lines of non-isolated singular points (or 
one-dimensional singular sets) L±. On the line L+ the singular points are either 
saddles or local sinks, and on the line L_ the singular points are local sources or 
saddles. The dynamics is very simple due to the existence of (two) monotonically 
increasing functions. A full analysis is given in |5B]. Henceforward, let us consider 
the case A > 0. 



4.2.1 Models with Positive Central Charge Deficit 

In the case where the central charge deficit is positive, A > 0, we choose the 

normalization 

/I \ ^/2 

e={3h^ + -p + N + 2A) . (4.17) 



2' 
The generalized Friedmann constraint equation ( [4.14| ) now takes the simple form 



^ = 1 (4.18) 



e2 



and may be employed to eliminate ip. Since by definition e > 0, specifying one of the 
roots ip/e = ±1 corresponds to choosing the sign of ip- However, it follows from the 
definition in Eq. ( [4.15|) that changing the sign of ip is related to a time reversal of the 
dynamics. In what follows, we shall consider the case ip/e = +1; the case ip/e = — 1 
is qualitatively similar. 



Introducing the new normalized variables 

/3h 

' '^- o 2' 

e 2e^ e 



/i = ^— , z/ = ^, X = - (4.19) 



and a new dynamical variable 



d 1 d 



dT y^e dt 
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(4.20) 



transforms Eqs. ([4.1CI|) - ([4.13|) to the three-dimensional autonomous system: 



dfi 
dT 
dv 

It 

dX 
'dT 



v + 






/i 



-2z/ 



/^ + 



U^ 



Vs 



fi 



Ts'^'-^' 



A) 



(4.21) 
(4.22) 
(4.23) 



The phase space variables are bounded, < {/x^, z/, A} < 1, and satisfy /i^ + z/ + A < 1. 
The sets z/ = and A = are invariant sets corresponding to p = and [3 = 0, 
respectively. In addition, /i^ + z/ + A = 1 is an invariant set corresponding to A = 0. 
We note that the right-hand side of Eq. ( |4.23| ) is positive-definite and this simplifies 
the dynamics considerably. 

l-( |4.23| ) consist of the isolated singular 
A = 0, (4.24) 



The singular points of the system 
point 

C : fi = u 

and the line of non-isolated singular points 

V:z/ = 0,A = 



fi^ [ji arbitrary). 



(4.25) 



The eigenvalues associated with C are Ai = I/a/S, A2 = 2/-\/3 and A3 = 0. Since C is 
an isolated singular point, it is therefore non-hyperbolic. However, a simple analysis 
shows that it is a global source. The eigenvalues associated with V are: 



Ai 



-2 /i + 



V3, 



A2 



2 
73 



(4.26) 



and the third eigenvalue is zero since ^ is a one-dimensional set. Therefore, on V the 
singular points are saddles for ji G [—1, — 1/v^) and local sinks for /i G (— 1/v^, 1]. 

It is also instructive to consider the dynamics on the boundary corresponding to 
A = 0, since the case iV = is of physical interest in its own right as a four-dimensional 
model. In this case the ODEs reduce to the two-dimensional system: 



dT ./5 V ^ 



dv 
dT 






'^Ts" 



The singular points and their corresponding eigenvalues are: 



C : 



fx 



V 



0; Ai 



1 



v^' 



A, 



(4.27) 
(4.28) 

(4.29) 
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S: ^ = -l,u = 0; X, = - — ,X, = 2\l--^j (4.30) 

A: /i = l,i. = 0; Ai = --|,A2 = -2('l + -^V (4.31) 

Point C is a non- hyperbolic singular point; however, by changing to polar coordinates 
we find that C is a repeller with an invariant ray 9 = tan~^(— -\/3). The saddle S and 
the attr actor A lie on the line V. 

4.2.2 Discussion 

Let us briefiy summarize the dynamics in the case A < 0, which was studied in 
p3| . When the modulus field is frozen the universe contracts from a singular initial 
state. The orbits are past-asymptotic to a dilaton-vacuum solution, and so the 
axion is negligible and the kinetic energy of the dilaton is dominant. As the collapse 
proceeds, however, the axion becomes dynamically more important and eventually 
induces a bounce. In the case of vanishing A, these equations imply that the future 
attractor would correspond to a dilaton-vacuum solution. However, the combined 
effect of the axion and the central charge deficit is to cause the universe to evolve 
towards a singularity with a —>■ +oo in a finite time. This behaviour is different from 
that found when no axion field is present in which there is no bounce []86| , |90| . The 



inclusion of a modulus field leads to a line of sources and a line of sinks. The axion 
field is dynamically negligible in the neighbourhood of the singular points. Moreover, 
a bouncing cosmology is no longer inevitable and there exist solutions that expand to 
infinity in a finite time. The solutions are asymptotic to the dilaton-moduli-vacuum 
solutions near the lines L±. 

In the case of a positive A discussed here, the isolated singular point C corresponds 



to the 'linear dilaton-vacuum' solution |^6|, |91|. When the modulus is frozen, all 
trajectories evolve away from C towards the point A and approach a superinfiationary 
dilaton-vacuum solution defined over t < 0. Some of the orbits evolving away from C 
represent contracting cosmologies and the effect of the axion is to reverse the collapse 
in all these cases. For the rolling modulus solutions, the orbits tend to dilaton- 
moduli-vacuum solutions as they approach the attractors (the sinks on V) . As in the 
case of a negative central charge deficit, the critical value /x^ = 1/3 corresponds to the 
case where d^ = /3^, representing the isotropic, ten-dimensional cosmology [a = (5) 
and its dual solution [a = —13). In the latter solution, the ten-dimensional dilaton 
field, $ = $ + 6/3, is constant. The other boundary of V is the point A representing 
the case where the kinetic energy of the modulus field vanishes. We note that the 
qualitative behaviour of models with ?/; < is similar. 
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4.3 Qualitative Analysis of the RR Sector 

When a cosmological constant is introduced into the matter sector of Eq. ( [4.4|) and 
the central charge deficit vanishes (i.e., A = 0), the FE are given by 

a = d0 + 0' - 3d' - 6/52 - ^ARe'^+3" (4.32) 

^ = 3^2 + 6/32 + Iarc'^+S" (4.33) 

a = -(vJ + 6d)(T (4.34) 

(3 = $ip (4.35) 

and the generahzed Friedmann constraint equation takes the form 

3d' - 0' + 6/3' + ^(T2e''^+'^" + ARe'^+2° = 0. (4.36) 
Eqs. ( [4.32D -( p:.36D may be simphfied by introducing the new time coordinate 

£^g~(^+3a)/2l (4.37) 

dt dt ^ ' 



and employing the generalized Friedmann constraint equation ( [4.36|) to eliminate the 
axion field. The remaining FE are then given by 

a" = y." - V + laV - 6/3'2 - ^Ar (4.38) 

^" = 3a" + 6/3" - ^^" - ^aV + ^Ar (4.39) 

P" = ^/3' (^' - 3a') , (4.40) 

where in this section a prime denotes differentiation with respect to t. 



4.3.1 Positive Cosmological Constant 

When Ar > 0, we express Eqs. ( [i.38| )-( |4.40| ) as an autonomous system by defining 



h = a\ ^ = if\ N = (3'. (4.41) 

Eq. ( [4.36| ) then implies that 

ifj^ > 3/1^ + 6N^ + Ar > (4.42) 
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and consequently we may normalize using ip. We therefore define 

^^ (4.43) 

(4.44) 

(4.45) 

(4.46) 

and assume that ip > [tp < is again related to time- reversal). 
The resulting three-dimensional autonomous system is: 

- \/3] (4.47) 



X = 


ip 


y = 


6N^ 


z = 


Ar 
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1 d 


de 


^£ 



dx 

de 

dy_ 

de 

dz 

de 



(x + V3)[l- 


2 T 1 r 

- X — y — z\ -\ — z[x 


2y([l-x^ 


-y-z] + -zj 


2zl[l-x'^ 


-y-z]--ii-z- 



(4.48) 
VSx) \ . (4.49) 



It follows from the definitions (|4.43| )- (|4.45| ) that the phase space is bounded with 



< {x'^,y, z} < 1, subject to the constraint 1 — x'^ — y — z>0. The invariant set 

1 — x"^— y — z = corresponds to a zero axion field. The dynamics of the system 
( [4.47] )- (|4.49| ) is determined primarily by the dynamics in the invariant sets y = 



and z = 0. These correspond to a zero modulus field and a zero Ar, respectively. 
The dynamics is also determined by the fact that the right-hand side of Eq. ( [4.48|) 
is positive-definite so that y is a monotonically increasing function. This guarantees 
that there are no closed or recurrent orbits in the three-dimensional phase space. 

4.3.2 Four-Dimensional Model 



In the invariant set y = 0, where the modulus field is trivial, the system ( [4.47|) -( p:.49|) 
reduces to the following plane system: 

dr 1 

_ = (x + V3)[l-x^-z] + -z[x-V3] (4.50) 

^ = 2z[[l-x'-z]-^-{l-z-V3x)Y (4.51) 

The singular points and their associated eigenvalues are given by 

Si: x = -l,z = 0; Ai = 2(^3-1), A2 = -(1 + Vs) (4.52) 

^2: x = l,z = 0; Ai = -2(^3 + 1), X^ = {V3 - I) (4.53) 



F: x = -^,z=if; Ai,2 = - ± -V23T. (4.54) 
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The points Si and 5*2 are saddles and F is a repelling focus. 

In the invariant set 1 — x"^ — z = 0, corresponding to the case of a zero axion field, 
Eqs. ( [4 .501) and ( |4.51| ) reduce to the single ODE 

§ = i (1 - x^) (x - v^) . (4.56) 

which can be integrated to yield an exact solution (in terms of G time). 

4.3.3 Ten— dimensional Model 



In the full system (|4.47|) - (|4.49|) with a non-trivial modulus field, there exists the 



isolated singular point F : x = —l/{3^/3),y = 0,z = 16/27 with the associated 
eigenvalues Ai 2 and A3 = 4/3. This implies that F is a global source. All of the 
remaining singular points belong to the one-dimensional set 

W:y=l — x'^,z = {x arbitrary) (4.56) 

and the associated eigenvalues are given by 

Ai = -2v^ fx + -^ j, A2 = v^U--^j, (A3 = 0). (4.57) 

W lies in the invariant set z = on the boundary y = 1 — x^. Points on W with 
X E (— 1/v^, 1/v^) ai'e local sinks, while the remaining points are saddles in the 
full three-dimensional phase space. (In the invariant set z = singular points with 
X E [—1, — 1/v^) are repelling and those with x E (— 1/\/3, 1] are attracting). 



We note that there exists an exact solution of Eqs. ( |4.471 )- (|4.49| ) with non-trivial 

modulus field with 

1 

X = 1= = constant, (4.58) 

3^3 ^ ^ 

and 

--f(^-i). («^) 

whence 

^ = - (- - -) . (4.60) 

de 4 V 277 ' ^ ' 

which can be integrated in terms of G time explicitly. 

Finally, we could consider the case Ar < 0. The generalized Friedmann constraint 
equation ( |4.36| ) implies that 

^^ - Ar > 3/i^ > 0, (4.61) 

and we could therefore normalize using y^xjj'^ — Ar. The resulting dynamical system 
can be analysed in a manner similar to that above (see p6|). 
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4.3.4 Discussion 

The mathematical structure of the dynamics is much richer when there is a positive 
RR charge. The dynamics of the system ( [4.50| )-( p:.51j ) is of interest from a mathemat- 
ical point of view due to the existence of the cyclic behaviour. The orbits are future 
asymptotic to a heterodinic cycle, consisting of the two saddle singular points Si and 
5*2 and the single (boundary) orbits in the invariant sets z = and 1 — x"^ — z = 
joining Si and 5*2. The former set corresponds to the zero Ar solution and the latter 
to the solution with constant axion field. In a given 'cycle' an orbit spends a long 
time close to 5*1 and then moves quickly to 5*2 shadowing the orbit in the invariant set 
z = 0. It is then again quasi-stationary and remains close to the singular point 5*2 be- 
fore quickly moving back to 5*1 shadowing the orbit in the invariant set l — x"^ — z = 0. 
We stress that the motion is not periodic, and on each successive cycle a given orbit 
spends more and more time in the neighbourhood of the singular points 5*1 and 5*2. 
The exact solution corresponding to the singular point F is power-law and is defined 
over the range — oo < t < 0, and represents a cosmology that collapses monotonically 
to zero volume at t = in which the curvature and coupling are both singular |]^ . 



The cyclical nature of the orbits can be physically understood by reinterpreting the 
axion field in terms of a membrane. Since the axion is constant on the surfaces of ho- 
mogeneity, the field strength of the two-form potential must be directly proportional 
to the volume form of the three-space. If the spatial topology of the universe is that 
of an isotropic three-torus, the axion field can be formally interpreted as a membrane 
wrapped around this torus |Q. All orbits begin at F. As the universe collapses, the 



membrane resists being squashed into a singular point. This forces the universe to 
bounce into an expansionary phase. The cosmological constant then dominates the 
axion field as the latter's energy density decreases. The subsequent effect of the cos- 
mological constant can be determined by viewing Eq. (f4.4|) in terms of a Brans-Dicke 
action, where the coupling parameter between the dilaton and graviton is given by 
5^; the late-time attractor then corresponds to a dilaton-vacuum solution. 



UJ 



In effect, therefore, the cosmological constant resists the expansion and ultimately 
causes the universe to recoUapse and asymptotically approach the saddle point 5*1. 
On the other hand, the collapse causes the axion field to become relevant once more 
and a further bounce ensues. The process is then repeated with the universe undergo- 
ing a series of bounces. The orbits move progressively closer towards the two saddles, 
5*1^2, and spend increasingly more time near to these points. This behaviour is related 
to the fact that the kinetic energy of the shifted dilaton field increases monotonically 
with time, since Eq. ( [4. 331 ) implies that (/3 > 0. It would be interesting to consider 
the implications of this behaviour for the pre-big bang inflationary scenario 



When a modulus field is included {y ^ 0), F still represents the only source 
in the system. The orbits follow cyclical trajectories in the neighbourhood of the 
invariant set y = and they spiral outwards monotonically around the orbit which 
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corresponds to the exact solution given by Eqs. (f4.58|) -( ^^60D , with dy/dQ > from 
Eq. ( [4 .481) . After a finite (but arbitrarily large) number of cycles the kinetic energy 
of the modulus field becomes more important until a critical point is reached where 
it dominates the axion and cosmo logical constant. The orbits then asymptote toward 
the dilaton-moduli-vacuum solutions corresponding to sources on the line W. 

A finite number of modulus fields and shear modes can be introduced into the 
model by defining the variable iV^ in Eqs. ( [4.41| ) and ( [4.44| ) via iV^ = J27=i ^h O^'" 



bits with non-trivial modulus field or shear term 'shadow' the orbits in the invariant 
set ?/ = at early times and undertake cycles between the saddles on the singular set 
W close to 5*1 and 5*2. These saddles may be interpreted as Kasner-like solutions 0; 
the orbits thus experience a finite number of cycles in which the solutions interpolate 
between different Kasner-like states. This is perhaps reminiscent of the mixmaster 



behaviour that occurs in the Bianchi type VIII and IX cosmologies R 97|. Mixmaster 



oscillations also occur in less general (i.e., lower-dimensional) Bianchi models with a 



magnetic field ||12[ or Yang-Mills fields [^. It is interesting to note in the string 



context that mixmaster behaviour also occurs in scalar-tensor theories of gravity in 
general and in Brans-Dicke theory in particular |^. Unlike the mixmaster oscilla- 



tions, however, the orbits eventually spiral away from y = 0, although there are orbits 
that experience a finite but arbitrarily large number of oscillations. 

Finally, we observe that all of the exact solutions corresponding to the singular 
points of the governing autonomous systems of ODEs are self-similar since in each case 
the scale factor is a power-law function of cosmic time. Therefore, exact self-similar 
solutions play an important role in determining the asymptotic behaviour in string 
cosmologies. However, we note that not all of the string solutions are asymptotically 
self-similar due to the existence of the heteroclinic cycle in the RR sector with trivial 
modulus field. 
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